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The activities during this reporting period (July 1 - December 

31, 1984) are described as follows. 

(1) Nonlinear wave drag 

The linear lifting-surface theory can predict reasonably 
well the lift in supersonic flow, but the drag is usually 
underpredicted. To remedy this deficiency, a simple method 
for calculating the nonlinear wave drag was developed. In 
the method, the calculated sectional drag is modified by 
adding the difference between the exact two-dimensional (2-D) 
and the linear 2-D wave drag at the calculated sectional 
lift coefficient. The results in most cases show good 
improvement in the supersonic drag prediction. All supersonic 
cases which were calculated previously were re-calculated. 

Two examples are presented in Figs. 1 and 2. 

(2) Restructuring of the VORCAM code 

At the suggestion of NASA users, the VORCAM code was 
modified for the Fortran 77 language. Its input stream 
was also rearranged. The work has been completed and has 
passed different tests so far. The only remaining task is 
the revision of user’s manual. 

(3) Boeing free-sheet code 

Great effort and time have been spent in adapting the 
Boeing code to the computer system at KU. All CDC special 
features in the code have been replaced with standard Fortran 
algorithms. This was accomplished with the assistance of some 
CSC personnel at the Langley Research Center. At this time, 
the code has passed successfully through compilation. Currently, 
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execution test is under way. 

(4) Investigation of asymmetric vortex separation 

To account for the asymmetric vortex separation, additional 
vortex distribution on the wing is needed to satisfy the boundary 
condition. Since the present preliminary investigation is 
based on the Smith's code of slender wing theory, this 
additional vortex distribution was introduced in his trans- 
formed plane (Z*) , instead of the physical plane (Z) . In 
the transformed plane, the wing is a vertical slit as shown 
in Fig. 3. In a symmetrical vortex configuration, this 
additional vortex distribution should vanish. 

Let V(z*) be the velocity induced on the transformed wing 
by the asymmetrical vortex system. To satisfy the boundary 
condition, a vortex distribution ("/) is introduced such that 
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Eq. (1) can be inverted to give 
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In the thin airfoil theory, C is always zero because of 
Kutta condition. In the present case, the leading edge is at 
y* = 0 and z* = 0. It follows that '/(O) = 0. Hence 
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To integrate Eq. (2), let 
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To evaluate the two integrals in Eq.(2), note that 
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Using Eqs. (9) and (10), Eq. (2) can be rewritten as 


tfz*)” —r 21 y=Jh= < u > 
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Using the above formulation, it is expected that if the 
asymmetrical vortex system can be predicted by an inviscid 
method, the additional vortex distribution in Fig. 3 would 
play an important role. So far, several ideas have been 
tested with the revised code. The vortex cores are always 
moved in a direction to minimize the force based on a combined 
gradient and conjugate gradient method. 

a) Movement of vortex cores. 

It is found that to have a better convergence 
characteristics one vortex core should be moved first 
than the other, instead of moving both at the same time. 

b) Step sizes 

Instead of having the same step size for ^ and 
Zy. , it is preferred to have separate step sizes in 
some cases. 

c) Updating of wing boundary condition 

Updating of wing boundary condition represents 
another level in the iterative method of solutions. It 
Is not apparent when and how often the wing boundary 
condition should be updated. In the revised code, option is 
provided to update the wing boundary condition either 
at the beginning only or as often as practical. 
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d) Convergence criteria 

The most difficult condition to be satisfied is the 
force free condition for the vortex and its cut. Instead 
of specifying a small error as a convergence criterion, 
it is preferred to start with a larger error and decrease 
it gradually. The results of any run, whether converged 
or not, are always kept on disk files for subsequent 
possible rerun. 

Some preliminary results are described below. Note that 
the starting asymmetrical vortex shapes are generated from the 
symmetrical shapes at different angles of attack. The results 
in Fig. 4 are for a & of 200 deg. (an angle through which a 
sheet turns around the core). It is seen that the resulting shapes, 
though not as asymmetrical as expected, are not symmetrical either. 
The configuration is further compared with the symmetrical case 
in Fig. 5. Another indicator of asymmetry is the vortex distri- 
bution (V). Fig. 6 shows the /-distribution. Note that is 
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At the suggestion of Dr. Dick Barnwell of NASA Langley, A 
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is increased to 450 deg. so that the vortex sheet will wrap 

around the core further. In this case, if the starting shapes 

on both sides are too much different, the force-free condition on 

the vortex sheets to determine the strengths is frequently violated 

and the calculation will just stop. Therefore, the starting 

shapes are generated corresponding to a = 8 and 6 respectively. 

Note that for o( = 25 deg., a = 7.134. The resulting shapes after 

-4 

reaching an average vortex force equal to £ = 10 are compared 
in Fig. 7 and the additional vortex distribution (/) is shown 
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in Fig. 8. Note that large y near z* - ±1 implies large tangenial 
velocity near the center of the physical wing. A tuft-grid 
survey of the flow field for this wing is presented in Fig. 9. 
Unfortunately, the survey was made at a station behind the wing 
trailing edge, instead of over the wing. 

A main conclusion from this study is that perhaps because of 
nonlinearity, the solution appears to be nonunique. One possible 
mechanism of vortex asymmetry has been "crowding" of two vortices. 
This is being investigated at the present time. 
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Figure 4 Vortex Shapes for 86.5-deg Delta Wing at a = 25 deg. 


Converged Symmetrical Shapes 

Final Shapes from Asymmetrical Configuration 
Shown in Figure 4 
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Figure A tuft-grid survey at about 0.2c r behind the trailing edge of a flat-plate delta wing with an 86.5° sweep of the leading edges for 

a = 25° and p = 0°. 
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